Perfect fluids and generic spacelike 
singularities 

Patrik Sandin AND Claes Uggla^'*' 

^Department of Physics, University of Karlstad, 
S-651 88 Karlstad, Sweden 

January 6, 2010 



Abstract 

We present the 1+3 Hubble-normalized conformal orthonormal frame approach to 
Einstein field equations, and speciafize it to a source that consists of perfect fiuids with 
general barotropic equations of state. We use this framework to give specific mathe- 
matical content to conjectures about generic spacelike singularities that were originally 
introduced by Belinskii, Khalatnikov, and Lifshitz. Assuming that the conjectures hold, 
we derive results about how the properties of fluids and generic spacelike singularities 
affect each other. 
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1 Introduction 

Although the singularity theorems say little about the nature of singularities, the very defini- 
tion of a singularity implies that there exists a variable scale — the affine parameter distance 
from/to the singularity of a causal inextendible geodesic that is used to define it, further- 
more, the one dynamical input that goes into the theorems, the Raychaudhuri equation for 
the expansion 9^ also implies a variable scale given by the expansion itself, since 9 has unit 
(time)~^ (or, equivalently, (length)"^, since we set the speed of light c to one). In this pa- 
per we study detailed asymptotic dynamical aspects of generic singularities, and this brings 
the expansion and the coupling of the Raychaudhuri equation to the remaining Einstein's 
equations into focus. We will locate the singularity in the past and we therefore refer to it 
as a 'cosmological' singularity. Since we study asymptotic temporal developments, we con- 
sider timelike reference congruences for which > in the vicinity of the singularity, where 
9 — >■ +CXD asymptotically, i.e., we are interested in 'crushing' singularities. Furthermore, due 
to the 'cosmological' context we will replace 9 with the Hubble variable H which is defined 
as _ff — i 6* (note that it is common in FRW cosmology to refer to H^^ as a characteristic 
time scale, also known as the Hubble radius when referred to as a length scale). 

The asymptotic blow up of H suggests that we should asymptotically 'factor out' H, 
and thereby the associated variable scale, toward the singularity, preferably so that the two 
following desirable features are incorporated into the formalism: 

(i) Preservation of causal structure, since it is reasonable to believe that there is a close 
connection between causal structure and the nature of singularities. 

(ii) Adaption to scale-invariance, since there are many known as well as conjectured links 
between scale- invariant, i.e., self-similar, solutions and asymptotic properties of many 
types of singularities. 

The natural way to accomplish this is by means of a conformal transformation (satisfies 
(i)) with a conformal factor that involves H (factoring out of H) so that the key variables 
are (conformally) scale- invariant, i.e., dimensionless, and thus adapted to the properties 
of self-similar solutions, since such solutions are scale-invariant (satisfies (ii)). Hence we 
use a conformally Huhhle-normalized scale- invariant formulation based on the conformal 
transformation 

G = H^g ^ g = H-'^G, (1) 

where we assume that H > in the vicinity of the singularity; g is the physical metric, which 
like H^^ naturally carries dimension (length)^, and hence it follows that the unphysical 
metric G is dimensionless. Because of this, scalars constructed from G take constant finite 
values for self-similar models that admit spacetime transitive homothetic symmetry groups. 
This leads to a major advantage: Asymptotically bounded variables for a system of coupled 
regularized field equations. 

We also find it advantageous to express the field equations as a system of first order 
partial differential equations. A natural way to do this within the conformally Hubble- 
normalized scale-invariant context is to use the Hubble-normalized Conformal OrthoNormal 
Frame approach (subsequently shortened to the acronym CONF). In this approach one 
chooses a frame field that is orthonormal to the dimensionless metric G, and not to the 
physical metric g, i.e., we introduce Hubble-normalized conformal orthonormal vector fields 
da that are dual to $7°, i.e., ($7°, db) = S°-b, such that 

g = H-''G = H-^Tjabn^^'' , (2) 

where rjab = diag(— 1, 1, 1, 1), and a,b — 0, 1, 2, 3. Furthermore, since we are interested in 
asymptotic temporal behavior, we let da be tangential to the reference congruence, i.e., 
do (X d/dx°, where is the time coordinate along the reference congruence (see Ap- 
pendix[A|. This naturally leads to the 1-1-3 Hubble-normalized CONF formulation, which is 



^In the case of timelike geodesies; in the null geodesic case an analogous equation plays a similar role. 
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a specialization of the 1+3 CONF formulation, introduced in [1 , to a conformal factor re- 
lated to H as described above; we present the 1+3 Hubble-normalized CONF field equations 
in Appendix El 

In this paper we consider a source that consists of several perfect fluids. The i:th perfect 
fluid yields a stress-energy tensor component, 

to the total stress-energy tensor, T"'' — ^ ■ T^^, where and are the energy density 
and pressure, respectively, in the rest frame of the «:th fluid, while M^j-j is its 4-velocity; 
throughout we assume that > 0. It is natural to make a 1+3 split of {t°j-j w.r.t. the 
vector held that is tangential to the reference congruence, and introduce a 3- velocity u"-^ 
according tcH 



«W =r(,)K + t;^,)); Uavf^^^O, r(,) = l/^l-i;2). (4) 

The i:th fluid is, apart from its 3-velocity, conveniently characterized by its energy-density 
w.r.t. u"', jO(j), which is defined in terms of /5(i) and w^^^ according to 

Pi^) =^{i)G+ Pi^), G^' = 1±W(,)W(,), w^i)^—. (5) 

P{i) 

Throughout, we are going to assume that the perfect fiuids satisfy barotropic equations of 
state, i.e., = P(i){p(i)), and hence it;(i) (/5(j)); special cases of interest are dust, w = 0, 
radiation, w = i, and stiff fluids, w = 1. 

To conform with standard convention in cosmology, we Hubble-normalize as follows 

where H is the Hubble variable associated with u'^. A 1+3 irreducible Hubble-normalized 
decomposition of the stress-energy tensor, see Appendix |3 yields 

= F(.)=u;(,)f7(,) + i(l-3«;(,))g««^,); H« = Q«««, (7) 

where Q^^j , P(i) , H^*^ are the Hubble-normalized components of the energy flux, pressure, 
and stress tensor, respectively (w.r.t the temporal reference congruence), and hence the 
Hubble-normalized stress-energy tensor of the i:th fluid is characterized by ^(i), v"^y and 
{a, (3 = 1,2,3). 

The 1+3 Hubble-normalized CONF formulation of the field equations for / fiuids involve 
the following quantities, see Appendix [K[ 

Hubble-normalized frame variables: {M, Ada, Ea^}- 

Hubble-normalized connection/commutator variables: {VF", If", i?", Sq.^, Aa, Najs}. 
Hubble-normalized perfect fluid variables: ^"i): ''^(i)j ■•••i ^(/)}- 

Apart from the already described fluid quantities, the quantities A4 and Ada are the Hubble- 
normalized threading lapse and shift functions respectively, while Ea^ are the Hubble- 
normalized spatial frame components. The quantities W"', U°', Sq/j describe the vorticity, 
acceleration, and shear of the Hubble conformal reference congruence, while R" describes 
the rotation of the spatial frame w.r.t. a Fermi frame in the space orthogonal to the refer- 
ence congruence. Finally Aa and Najs gives the commutator functions (or, equivalently, the 
spatial Hubble-conformal connection coefficients) of the Hubble-normalized spatial frame. 



■^One reason for why this is convenient is that since the components of the 3-velocity in the orthonor- 
mal frame of g are dimensionless they coincide with the 3-velocity components of the conformal 4-velocity 
in the Hubble-normalized frame of G. 
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The outline of the paper is as fohows. In the next section we formulate two conjectures, 
originally introduced by Belinskii, Khalatnikov, and Lifshitz, in terms of the above Hubble- 
normalized variables. In Section [3] we explore the consequences of these conjectures and 
derive a number of results concerning the past stability and instability on the so-called 
silent boundary; in particular we present the past attractors on the silent boundary (note 
that since there exists a dynamical one-to-one correspondence between the silent boundary 
and the spatially homogeneous models, it follows that our results also pertain to the latter 
case) . Then the stability of these results are investigated in the context of the full state space 
in Section m where we also discuss possible temporal gauge choices, and in particular if it 
is possible to use fluid congruences as temporal reference congruences to describe so-called 
asymptotically silent and local singularities; we find that this is only possible if there exist 
fluids with a sound speed that is equal to or larger than the speed of light. We conclude with 
a summary and some remarks about our results in Section [5l together with some comments 
about some open issues. Appendix [X] establishes conventions and notation by giving the 
l-|-3 Hubble-normalized CONF field equations explicitly. Finally we describe a number of 
important subsets in Appendix |BJ 



2 BKL conjectures 

In [2] p656 Belinskii, Khalatnikov, and Lifshitz made the following important conjecture: 

Conjecture, "...in the asymptotic vicinity of the singular point the Einstein equations are 
effectively reduced to a system of ordinary differential equations with respect to time: the 
spatial derivatives enter these equations 'passively ' without influencing the character of the 
solution. " 

In the present 1+3 Hubble-normalized CONF framework we reformulate this conjecture 
in terms of two intertwined conditions: 

Conjecture (1). 

(a) Asymptotic surface formation : lim {Aia, Wa, Ua, ^a) = 0, 

< Ci < lim M<C2<oo. (8a) 
(h) Asymptotic locality condition : lim {E^'^^ 5a X) — 0, (8b) 

x'-^-^ — oo 

where Ci,2 = const and X = (M, Ma, Wa, U", R" , Sa/3, A^, N^fi, rj(i), v'^^y w^^^), 
where I denotes the number of perfect fluids. 

Condition (a) implies that the spatial frame is asymptotically hypersurface forming; 
in addition Ua — Q implies that the timelike congruence is conformally geodesic, which 
amounts to an inverse mean curvature flow for the original physical spacetime, while r^ — Q 
implies that a foliation is a constant mean curvature foliation in the physical spacetime. 
Furthermore, Ua — r^ — Q implies that the reference timelines are geodesies in the original 
physical spacetime. Finally, the condition that M is asymptotically bounded implies that 
the reference congruence asymptotically gives rise to a foliation that yields a simultaneous 
bang function when — cx)|j 

The field equations for the 1-1-3 Hubble-normalized CONF- variables for a source that 
consists of several perfect fluids are given by equations ([73|) . ([75|) . and ([79| in Appendix lAl 
This system admits an invariant subspace called the silent boundary., see Appendix |B1 which 
is characterized by 

(A^„, W^a, r„, =0. (9) 

^This follows from that we can reparameterize the congruence according to = exp(x''), so that the past 
singularity occurs at = 0. However, we believe that the condition on M can be weakened, which, however, 
we have refrained from doing in order to keep the discussion reasonably simple. The central restriction is to 
choose a reference congruence that asymptotically yields a spacelike foliation that has a simultaneous bang 
function, cf. the discussion about the synchronous gauge in |2]. 
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As discussed in Appendix 1X1 the equations on this subset form a coupled system of ODE that 
is identical to the system that describes the dynamics of spatially homogeneous Bianchi mod- 
els. If conjecture (1) holds along a timeline we may say that the dynamics become asymp- 
totically local, since the dynamics then is asymptotically described by the silent boundary 
on which the dynamics for a timeline is governed by what happens along the timeline alone, 
and hence it would perhaps be more appropriate to refer to the silent boundary as the local 
boundary^ Finally, a singularity that obeys 'the locality conjecture' (1) will be referred to 
as an asymptotically local singularity. 

The dynamical relevance of the silent/local boundary depends on if the conditions in 
conjecture (1) holds. From now on we will assume that this is the case and derive the 
consequences of this assumption. A necessary condition for the dynamics of a timeline to 
approach the silent/local boundary is that i?Q,* — >■ toward the past singularity, which is 
equivalent to that the conformally Hubble-normalized contravariant spatial 3-metric 

SQrj ^ gafS (^Q^ 

tends to zero. Due to (|73c|) . ■t'*-' satisfies the equation 

do^G'^ = 2{qS"^ - J:"'^)E^' EpK (11) 

The vanishing of Ea to the past is equivalent to the condition that the time integral to 
the initial singularity of the eigenvalues of the matrix {qSa^ — Sq,^) negatively diverge for 
all of the eigenvalues, 

" — CO 

eig{q6j - T.J)Mdx° = -oo, (12) 

where describes some reference point along the timeline. It has implicitly been shown 
that the condition (IT2t is fulfilled for vacuum and orthogonal fiuid Bianchi models of type 
IX [HEllT], but this, of course, does not imply that it is true in general. However, we will 
in the following assume that condition (fT2|) holds and work out the consequences of that 
assumption, which leads to a consistent picture. Our assumption that there exist generic 
dynamics that is asymptotically described by the silent boundary suggests an analysis in 
two steps: 

1. Identification of the past attractoiH on the silent boundary. 

2. Perturbation of the past attractor in the full inhomogeneous state space to establish 
if it is stable or not. 

A proof that identifies the attractor and shows its stability in the full infinite dimensional 
state space amounts to a proof of a singularity theorem that concerns the details of a generic 
singularity. This is likely to be an extremely hard problem, and we will therefore only provide 
proofs about some aspects in the context of that our 'BKL-like' assumptions hold. 

There exists a second conjecture proposed by Belinskii, Khalatnikov, and Lifshitz [51 [TU] 
that is relevant in this context; the asymptotic 'matter does not matter' conjecture. 

Conjecture. For a typical cosmological model, the matter content is not dynamically sig- 
nificant near the initial singularity. 

In our case of a source of several perfect fluids we formulate this conjecture in terms of 
our variables as: 



■* Asymptotic silence is defined as the formation of particle horizons that shrink to zero size in all directions 
along any timeline that is not asymptotically null toward the singularity, thus asymptotically prohibiting 
communication. When the nomenclature silent boundary was introduced it was believed that asymptotic 
silence implied asymptotical local dynamics, however, the discovery of recurring spike formation [3l|4] shows 
that this is not the case. Although we expect that conjecture (1) holds for most timelines for an open set 
of solutions, the existence of recurring spike formation suggests that there exist generic singularities with 
special timelines with 'non-BKL' behavior for which conjecture (1) does not hold. 

^The past attractor of a dynamical system given on a state space X is defined as the smallest closed 
invariant set C X such that the a-limits of all p S X, apart from a set of measure zero, satisfy 

a{p) C A- E. 



3 PAST STABILITY AND INSTABILITY ON THE SILENT BOUNDARY 



6 



Conjecture (2). 

lim ritot = 0, (13) 
where the total Hubble-normahzed energy- density of the source is given by 

i 

Note that fitot — implies r2(i) = OVi, since > 0, and from this it follows that the entire 
total Hubble-normalized stress-energy tensor is zero. Hence 'the matter does not matter' 
conjecture (2) asserts that the Hubble-normalized stress-energy tensor asymptotically ap- 
proaches zero toward the singularity. As pointed out by BKL themselves in the context of 
one fluid, this is not to be expected for all equations of state, e.g., not for a stiff fluid. 



3 Past stability and instability on the silent boundary 

To proceed with step 1 we first give the equations on the silent boundary, which are obtained 
by restricting the full system of equations, (|73l) . fTS]) . and (f79|) in Appendix lAl to the silent 
boundary invariant subspace, {Ma, Wa, Ua, ra, E^) = 0. As discussed in Appendix [Bl 
the equations on the silent boundary are the same as in the spatially homogeneous case, and 
hence all results in this section also pertain to these models. 

3.1 Equations on the silent boundary 

Instead of the peculiar 3-velocity we find it useful to introduce w > and the unit vector 
c" — /v as variables. This leads to the following state vector on the silent boundary: 

S = iV"^) ® (17(1), ^^(1), cfi)) ® ... ® {^(i), «(/), cf,)). (15) 

Note that we have not included Ra in the state vector since there exists no evolution equa- 
tion for which is due to that Ra represents the freedom to rotate the spatial frame 
(nor M, which represents the freedom to reparameterize the reference timelines). On the 
silent boundary we have the following evolution equations and constraints that govern the 
dynamics of S. 

Evolution equations: 

do ^a0 = -(2 - q)^ai5 + 2e^\a ^p)s ~ ^Sap + 3H„/3, (16a) 

doAa^Fj^Asi, (16b) 

do iV"'^ = (3g(5^'" - 2i^^(")iV'3)'^, (16c) 

do 17 = (2g - 1 - 3w) + [(3u; -~\)va- ^afi + 2Aa] Q", (16d) 

dov = GZ^ (1 - v^) [3cl ~l-2clA^ cpv -Y^af} c°'cl^] v, (16e) 

doCa = -[Sa'^ - CaC^ji^p'' + V Ap + tfs^^ {Rg + W iV/ C,) C^]. (16f) 

Constraint equations: 

= 1 - - Ofc - 17, (17a) 

= {SSa'' A0 + eas'' N^p) - 3Qa, (17b) 

= ApN^a, (17c) 

where 

Fa^ = q6a^ -Ya^ -ea^jiW^ + R^), g = 22^ + i (17 + 3F) , (18a) 

35a0 = B^al3) + 2e'\a Np)S A^ , Ba^ = 2Na^ N^ - N"' ^ Nap, (18b) 

37^= -iS"„-6A2, f!^, = _i37^, (18c) 

G- = l-c^^ c2 = |, (18d) 
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where = ^S„^S"^, and where can be interpreted as the speed of sound when non- 
negative. Note that it is the complete stress-energy-momentum objects that appear in (|16a|) . 
([T7a|) . (I17b|) . and ([T5a|) . while the perfect fluid equations (|16dl) and (fTBe)) . describe the 
dynamics of an individual perfect fluid component, where we have dropped the index (i) 
to avoid cluttered notation. To obtain the perfect fluid equations we have assumed that 
the Hubble-normalized interactions between the different fluids are asymptotically zero, see 
Appendix [A) It follows from (|16e|) that u = is an invariant subset and so is u = 1 when 

^ 1, i.e., when the equation of state of the fluid is not stiff. Remarkably, w does not 
appear in the peculiar velocity equations (|16cp and (jl6f|) . nor do 17 and q — the equation of 
state enters via only, and thus a general barotropic equation of state leads to formally the 
same expressions as that of a linear equation of state! However, in general is a function 
of a suitable matter variable, e.g. c^(p), while — w — const in the linear case. Moreover, 
the equation for the peculiar velocity direction Cq., i.e. (jl6fj) . contains neither nor w, i.e., 
it contains no direct coupling to the equation of state at all! 

For completeness we here give the evolution equation for the peculiar velocity Va on the 
silent boundary: 

dovo. = GZ' [(1 - v^){3cl -l~clAl' vp) + (1 - cl){Al' + E/ t;^) vp] 

- + eJ-< [R^ + vs)] vp -A^v\ (19) 

It is of interest to also give the evolution equations for p and p, for a fluid component, on 
the silent boundary (i.e., let {Ea,Ma, Wa, Ua, '''a) = in the equations for these objects): 

do{lnp) = -{l + w)G-^[3-2A^v°' + {v^ + J:^^v"'v^)], (20a) 
do (Inp) - -(1 + w)GZ^[3 - 2A„t;" - {v^ + v"v^)], (20b) 

where we again have dropped the index (i) in (fTO|) and ([201 to avoid cluttered notation. 
3.2 Past evolution on the silent boundary 

Proposition 3.1. The past asymptotic limit in a regime where H > resides on the type 
I — VII part of the silent boundary if the strong energy condition is asymptotically fulfilled. 

Proof. We have 

aodet(A^„^) = 3gdet(A^a^), (21) 

on the silent boundary, where 

g = 2E2+ i(r!t,t+3Ptot). (22) 

If the strong energy condition J^tot + 3Ptot > holdsj^ then q > 0, and q = only when 
ntot + 3Ftot = and = 0, but then 

dl det(Ar„^)|,=o = 0, dl det(iV,^)|,=o - 2 [""S^^ %s] det{N^p), (23) 

where ^S''^ ^S-ys > when det{Nai3) 7^ 0; it follows that 

det(Af„^) ^ (24) 

toward the past singularity. Thus the past asymptotic limit of the dynamics must reside on 
the det(A'Q,^) — subset, i.e., the Bianchi type I- VII part of the silent boundary. □ 

Corollary 3.1. 

lim |E„/3|<2. 

®It is likely that there exist generic solutions with only a positive cosmological constant as source (with 
S^tot + 3Ptot = —2Qtot < 0) that asymptotically behaves as generic vacuum solutions with asymptotically 
silent and local past singularities, and hence it should be possible to relax the condition Qtot + 3Ptot > 0, 
but for simplicity we refrain from doing this. 
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Proof. det{Nai3) — implies that ilk > 0, which, together with the Gauss constraint 1 — = 
flk + ft > 0, yields 

< 1 ^ -2 < < 2. (25) 

□ 

On the silent boundary 

doA^^2{q6j ~J:J)A"Af,, (26) 



and hence, assuming the validity of 'the locality conjecture' (1), and thereby that Eq. (|12l) 
holds, 

A„ -j> (27) 

toward the singularity, i.e., the past attractor has to reside on the subset that consists of 
the union of the class A {Aa — 0) type I, II, VIq, and VIIq subsets on the silent boundary. 
On the class A part of the silent boundary (|16eP reduces to 



dov ^ GZ^ (1 - v^) {3cl c^'c^) V. (28) 



Corollary (|3.1|) and Ec^. (f28|) indicate that there is a bifurcation in the dynamics of the 
particular velocities of the fluids when = 1. We will therefore below distinguish between 
three main cases, based on the asymptotic properties of the equations of state: 

(i) There exists at least one fluid with an asymptotically ultra-stiff equation of state, i.e., 
Cg > 1, ui > 1 when — ?> — oo. 

(ii) All perfect fluids have asymptotic equations of state such that < 1, w < 1 when 
x'^ — s> — oo, except for at least one fluid which has an asymptotically stiff equation of 
state, i.e., = l,w = 1 when — oo. 

(iii) All perfect fluids have asymptotic equations of state such that < 1, w < 1 when — >■ 
— oo, i.e., all equations of state are softer than a stiff equation of state asymptotically. 

We will denote the three cases as the (asymptotically) ultra- stiff , stiff , and soft cases, re- 
spectively; as we will see^thcir past dynamics is associated with an increasingly complicated 
and challenging analysis |3 

To proceed we prove, under assumption (|12p . the following lemma: 

Lemma 3.2. 

/ eig{25j - i:J)Mdx° = -oo. (29) 

Proof. Case (i): As shown next, Eq,^ — in case (i), and hence the integral diverges. Case 
(ii) and (iii): Eq. ^ yields iltot > Ptot, which leads to the inequality q = 2Y? -t- \ {^tot + 
iPtot) = 2 - |(rJtot - Ptot) ~ ^flk < 2. This combined with Eq. (dH gives 

r — OO r — OO 

/ eig{2Sj-Ej)Mdx"= / eig[{2 - q)S J + {q6 J -J: J)] Mdx° < 

e\g{q5j - J:J)Mdx° = -oo. (30) 

□ 

Proposition 3.2. The past asymptotic state in case (i) is characterized by 

f^ultra-stiff — > 1; {^af!: Naj3, Wultra-stifl , ^{i)) 0, Vi ^ ultra-stiff. (31) 



^The physical status of an ultra-stiff equation of state can be questioned since Cg is larger than the speed 
of light, however, it is of interest for structural stability reasons to study sources with fluids with such an 
equation of state, moreover, in [111 , and references therein, the study of problems associated with ultra-stiff 
equations of state is motivated by considering broader theoretical contexts than general relativity. 
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Proof. Corollary p.ip and Eq. ([28|) give that \im^o^_aa Wuitra-stiff — 0. On the class A 
t^uitia-stiff = boundary Eq. (|16d|) yields 

5o In(r2(i)/r2ultra-stiff) = 3 (Wultra-stiff — 1) 

+ (G+)(:i[3(l - w^,)){l - vl)) + (1 + w^,^){2S^p - ^c.0X,)V^^^^] > 0, (32) 

where the (i):th fluid has a comparably asymptotic soft equation of state. For simplicity 
we have assumed that the ultrastiff fluid obeys an asymptotically linear ultra-stiff equation 
of state such that Wuitra-stiS = limpuitra-=titt-*-oo(w); in the case of several fluids with the 
same asymptotic ultra-stiff asymptotic equation of state, fiuitra-stiff, represent their total 
contributions. Since the r.h.s. of (1321) is strictly positive it follows that ri(i)/fiuitra-stiff ~> 
toward the past, and since iluitra-stiff is bounded, because of the Gauss constraint 1 — — 
f2fc — fitot — and the non-negativity of the energy densities and flk, this leads to that the 
ultra-stiff fluid(s) dominates toward the singularity, and hence r2(i) — >■ 0; thus the attractor in 
the ultra-stiff case (i) resides on the class A Bianchi type I - VIIq part of the silent boundary 
with Wuitra-stiff = 0, ft(^i) = 0, for all i except for the i associated with the ultra-stiff fluid(s), 
subset. This leads to that (If 6dp asymptotically yields 

^O^^ultra-stiff = — [3(Wultra-stiff ^ 1)(1 ^ ^^ultra-stiff ) + 4r2fe] f^ultra-stiff , (33) 

and hence, due to that iluitra-stiff < 1, asymptotically riuitra-stiff = 1 and ftk — 0, and 
thus, because of the Gauss constraint, — 0. That fifc = and = yield that the past 
attractor in the ultra-stiff case must reside on the isotropic type I subset or the isotropic type 
VIIq subset; in the latter case we can choose a Fermi frame in which Nap = diag(0, A^, iV), 
or cycle, which yields do N = qN = i(f + 3'u;uitra-stiff)A'', and hence N ^ 0, i.e., the past 
attractor is located on the isotropic type I subset, which is a frame independent statement; 
we will refer to the silent isotropic type I subset as the silent Friedmann subset J- . 

The above arguments are easily generalized to the situation when the most ultra-stiff 
equation(s) of state does not have a limit, but a lower bound > f ; one still obtains 

that the past attractor resides on T with fiuitra-stiff = 1, ^{i) = 0, i^uitra-stiff = 
even though q has no limit. □ 

Proposition 3.3. The past asymptotic state in case (ii) is characterized by 

q^2, {Na0, Wstiff, ^(t)) -> 0, Vi ^ stiff. (34) 

Proof. The analysis of the stiff case (ii) proceeds with similar arguments as in the proof of 
case (i) , but with the extra condition of lemma 13.21 This leads to that Wstiff — asymp- 
totically, and that =0 asymptotically for all fluids with equations of state that are 
asymptotically softer than the asymptotically stiff fluid(s). Hence the past asymptotic state 
resides on the imion of the class A Bianchi type I, II, VIq, VIIq subsets for a single orthogo- 
nal stiff fluid, where c"jj and V(^i) act as test fields, i.e., fields that do not affect the spacetime 
geometry but are affected by it. The past asymptotic dynamics for the single orthogonal 
stiff fluid case in Bianchi types I, II, VIq, VIIq is well known [UIT^, and from this it follows 
that the past attractor resides on the type I subset where flgtiff = ^^stiff, 9 = 2, where we 
have introduced the convention of using hats on purely spatially dependent, i.e., temporally 
constant, quantities. It therefore follows that the past attractor in case (ii) resides on the 
type I subset where 

^tot = f^stiff = f^stiff, ^{i) = 0, Ustiff = 0, q = 2; (35) 

we will refer to this subset as the silent Jacobs subset J (the exact solutions for a single stiff 
perfect fluid in Bianchi type I were first found by Jacobs CH 

We now turn to the behavior of ri(i) in the soft case (iii). For this case we have no proof, 
but we expect that the 'matter does not matter' conjecture (2) holds, and that ilj^j) — ^ for 
all i toward the past singularity, and that the past attractor hence resides on the vacuum 
subset r^tot — 0. The reason for the expectation that iltot = asymptotically is that 
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there exists evidence for that this happens when one has one fluid with a soft equation 
of state [5l |6l [14], and it seems reasonable that one can apply this result for each fluid 
individually; furthermore, in the vacuum case there exists evidence that the past attractor 
resides on the union of the silent vacuum type I subset, known as the silent Kasner subset 
K,, and the silent vacuum type II subset [3l [Sj [6l [TH [15] . Moreover, in two previous studies 
of tilted multi-fluid models of Bianchi type I [1^1 [TT] we presented evidence that indicated 
that the past attractor of the Bianchi type I models with two soft fluids resided on /C, and 
since we expect that JC plays a 'dominant' role in the asymptotic dynamics this gives further 
support for the claim that fltot — 0. 

From Proposition 13.21 in case (i). Proposition 13.31 in case (ii), and the 'matter does not 
matter' conjecture (2) in case (iii), it follows that asymptotically toward the past 

Qrot=0 and K,{ = 0, (36) 

in all cases, since ilj-jj — 0, for all i, except for the asymptotically 'dominant' ultra-stiff 
fluid(s) in case (i) and the asymptotically stiff fluid(s) in case (ii), but in those cases 
I'uitra-stiff = and VstiS = 0, respectively. 

For all Class A models with Q^^^ = and Utot = it is possible to simultaneously diago- 
nalize Nap and T,ap in a Fermi frame. The reason for this is as follows: In class A Qf^^ = 
leads to that the Codazzi constraint ()17bl) takes the form £as'* p S^^ = 0, which implies 
that Naj3 and are simultaneously diagonalizable for a given arbitrary value of {Nap 
transforms as a tensor density on the silent boundary under spatial frame rotations). Fur- 
thermore, the preservation of the simultaneous diagonalization during evolution is possible 
because n"j^ = 0, but it also requires that one uses a Fermi frame. We hence expect that it 
is possible to asymptotically diagonalize Sq^ and N^p in a frame that is asymptotically a 
Fermi framed so that 

Ra^O, =diag(Si,S2,E3), Nap = disig{Ni,N2,N3); Si -f E2 -f E3 = 0. (37) 

In all fluid cases, the silent Bianchi type I subset plays a prominent role, indeed, according 
to the previous analysis the past attractors for the ultra-stiff and stiff cases reside there, and 
we therefore now turn to this subset in more detail. 



3.3 The silent Bianchi type I subset 

As follows from the previous subsection for cases (i) and (ii), and as conjectured for case 
(iii), sources that consist of multiple perfect fluids lead to that the past asymptotic subset 
for Bianchi type I resides on the subset with r2(j) — OVi, except for the 'dominant' matter 

component (s) Jluitra-stiff = 1 (where also ^'ultra-stiff = 0) and Qstis = f^stiff (where also 
w^jff = 0) in cases (i) and (ii), respectively. This implies that the past asymptotic dynamics 
on the silent type I boundary resides on J ^ and /C, for the ultra-stiff, stiff, and soft cases, 
respectively, where K, constitutes the boundary of J in the stiff case. In S^^-space the fleld 
equations for Bianchi type I immediately lead to that these subsets are characterized by the 
following eigenvalues for Sq,^: 



T : 




= Eq. 


= 0, V a. 


E2 


= ^ 


i^tot — 1. 


J: 


Eq, 


= 


V a, 


E2 


= E2 = 


1 ^ ^^tot = 1 ^ f^stifl 


K. : 


Eq; 


= Eq., 


V a, 


e2 


= Y? = 


1 ^ Q.tot = 0. 



(38a) 
(38b) 
(38c) 



The eigenvalues Ec 
defined according to 



Eq, can be expressed in terms of the shape parameters pa , see |18| 



(El, E2, E3) = (3pi - 1, 3p2 ~ 1, 3p3 - 1) 



Pi +P2+ P3 



(39) 



*At least generically, the results in |15| show that some degrees of freedom only seem to be statistically 
suppressed, and that there may be a few timelines with different behavior; hence there may be some timelines 
with different asymptotic dynamics than that we presently describe. 
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where we have omitted the hats on the spatially dependent pa to conform with standard 
notation. In terms of the shape parameters, the past asymptotic states on the silent Bianchi 
type I subset for the three subsets are described by: 



J 



{Pi,P2,Pz) = i (1,1,1). (40a) 
pI+pI+pI = l-|f^stiff< 1. (40b) 
pI+pI+pI = 1- (40c) 



Even though flj^) = Vi (with the exception of r^uitra-stiff = 1 and r^stiff = f^stiff in cases 
(i) and (ii), respectively), the subsets J^, J, and /C also involve the equations for c".^ and 
which act as test fields, i.e., fields that do not affect the spacetime geometry but are 
affected by it, and hence a complete past asymptotic description also involves the asymptotic 
determination of these fields. For this purpose we use a shear diagonalized Fermi frame so 
that = diag(3pi — 1, 3p2 — 1, 3p3 — 1) and Ra = 0, and insert Aa = 0, Na/s — 0, which 
characterizes Bianchi type I, into the equations (|16f|) for Cq,; this leads to: 



doci = 3[(p2 -Pi)cl + {p3 -pi)cl]ci, (41a) 



^3J 

doC2 = 3[(p3 ~P2)cl + {pi -P2)cl]c2, (41b) 

^0 C3 = 3[(pi - P3)cl + {p2 - P3)cl] C3, (41c) 

where we again for simplicity have dropped the index (i). These equations, which decouple 
from the equation for can be treated as a separate dynamical system that satisfies the 
constraint CqC" = 1, i.e., we have a dynamical system on a sphere with unit radius, parame- 
terized by pi, p2, and p^. We note that this system is the same as that for v" when = 1, 
i.e., the dynamics for Ca is the same as for the extreme tilt subset v'^ — 1, which in [14 was 
examined by means of spherical coordinates in the case pi < p2 < Ps- 

Proposition 3.4. The past asymptotic state of the system (|4ip is given by: 

Case (i): Ca — 

Cases (ii) and (Hi): Let {a/3"f) = (123), or a permutation thereof, (a) If Pa < P/3 < P7, 
then Ca, C[3 0, — >■ ±1 when c^ ^ 0. (b) If Pa < Pp = P-y, then Ca — >■ 0, cp, c-y — >■ 

C/3 , C'y , + — 1 . 

Proof. In case (i) where (pi,P2,P3) — ^(l, 1, 1), it follows directly that d^Ca ~ and hence 

Cases (ii) and (iii) can be treated collectively. We first note that if <Pi3 < P7, where 
(a/37) — (123), or a permutation thereof, then the system (|4ip admits the invariant subsets 
C12 on which C3 = 0, and cycle, leading to a division of the sphere into six disjoint subsets 
with the subset C12, C23, C31 as boundaries, furthermore, the intersections of these subsets 
yield the fix points for which Ca = ±1, = = 0, (a/37) — (123), and cycle. If 
Pa = Pp ^ P7, where (a/37) = (123), and cycle, then the system (|4T|) also admits subsets 
when one of the components Ci, C2, or C3 is zero, but the subset Cap on which = reduces 
to a circle of fix points with Ca — Ca, cp = cp, + c^^ — 1, which we denote by C^^. 

If Pa < Pp < Pj, then doc^ > 0, and hence Ca 0. This reduces the system to the 
subset Cp~f where pp < p^ gives 5o c| > 0, d^c^ = f{cp)c^, where f{cp) < 0, and hence 
cp ^ and c-y ±1 toward the past when ^ 0. If Pa < Pp — P-y, then Ca — still holds, 
but in this case the system reduces to the circle of fix points C^^. □ 

Corollary 3.3. The past asymptotic peculiar velocity direction(s) c" of the test fields co- 
incide with the asymptotic eigenvector(s) of Yiap associated with the eigenvalue(s) Smax — 
max(Ei, S2, S3) — max(3pi - 1, 3p2 - 1, 3p3 — 1). 

Proof. This follows immediately from the proof of p.4p . □ 
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We now turn to the past asymptotic behavior for the pecuhar test speeds v on the type 
I subset. By regarding Cq. as time-dependent coefficients in the evolution equation for v, we 
can apply a theorem by Strauss and Yorke ^19 that implies that v is past asymptotically 
determined by the past asymptotics of c^. CoroUarv 13.31 then reduces Eq. (|28l) for v to 

dov = 3G;^(1 - v^){cl -_Pmax)w, whcrc Pmax = max(pi , P2 , Ps) ■ (42) 

Consequently v is monotonically decreasing (increasing) toward the past if > Pmax = 

1(1 + Smax) (C^ < Pmax = |(1 + ^max)), and hcnCC V ^ {v = 1), whilc w = -u if = Pinax, 

asymptotically toward the past; in these formulas refers to the asymptotic limit of cj. 
when /5 — > oo (for simplicity we assume that cj. has such a limit, however, many of our 
results are easily generalized to the case when has asymptotic bounds, but no limit). 

Case (i): Pmax ~ h and hence v ^ when > ^; v ^ v when = i; u — >■ 1 when 

3. 

Case (u): Eqs. 1^ and pOb)) yield that i(l + t) < p^ax < 5(1 + 2S) < 1, where 

E = Vl — Ostiff < 1. Hence < i ^ w — 1 toward the past; if > i there exist 
some Pmax values on for which w — > 0, some for which «—>?), and some for which 
V ^ 1, depending on if > Pmax, = Pmax, or cj. < Pmax (the smallest possible Pmax 
value is | and occurs when = 0). 

Case (iii): Eqs. ^ and (I40bp yield that | < Pmax < lEI Hence < | ^ w ^ 1 
toward the past. If > | there exist some points on /C for which v ^ and some for 
which w — > 1. 

It is of interest to note that 1 w — ?> everywhere on J and /cF°l 



3.4 Stability and instability of the type I subset on the silent bound- 
ary 

For the ultra-stiff case (i) it is easily seen that J" is a stable subset w.r.t. perturbations of 
Ea^, Aa, Nai3, Eq^, and hence there exists a past attractor in the full state space that 
resides on in this case, a statement that is also supported by the analysis in [11]. We 
therefore turn to the past attractor for the stiff and soft cases (ii) and (iii), respectively. 

To identify the past attractor subset on the silent boundary we next linearly perturb J 
and JC by using a Fermi frame in which the perturbed Bianchi type I subsets are expressed in 
a Fermi frame with diagonalized shear Eq,^ = diag(Ei = 3pi — 1, E2 = 3p2 — 1, E3 = — 1): 

A-^doAa\jx='2-^a=Hl-Pc), (43a) 

iV-i do Na\jx - 2(1 + E„) = 6pa where N^ = N^a, (43b) 

N-^ do Nap\jx = 2 - E^ = 3(1 - p^) where (a/37) = (123) and cycle, (43c) 

n-^do \jx - 3G;1 [(1 -w)il- «2) + {l+w)il- p„ax) v^] , (43d) 

where the above equations refer to separate components. The notation \ j^ic indicates eval- 
uation at J7 in the stiff case (ii), and at K, in the soft case (iii). Note that we have obtained 
the same form for the equations in the stiff case (ii) and the soft case (iii), since q = 2 in both 
cases. Moreover, in the case of rt(i) we have in addition inserted the type I past attractor 
value associated with J, K. for E^^ c" according to Corollarv 13.31 

Notably there are no equations for Eq,^ in (|43l) . The reason for this is that the stability 
analysis of E^^ depends on the choice of spatial frame. However, Eqs. (|43p hold for any 
spatial frame that admits E^^ = diag(Ei, E2, E3) and = as an invariant subset on 



®We exclude that pmax = 1; this is intimately connected with that the assumption II12I I holds. Note that 
Pmax = 1 has been proved to be excluded for the non-LRS Bianchi types VIII and IX cases (5] [G] I20| . 
'^'^ Except for the excluded points with Pmax = 1- 
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i7 and IC; furthermore, when projecting out the pecuhar velocities, which we wiU do in the 
reminder of this subsection, these sets form sets of fix points for the projected system of 
equations. Nevertheless, in the full state space, as well as on the silent boundary where in 
general Aa ^ 0,Qa O,!!^^ ^ 0, the shear cannot be diagonalized in a Fermi frame, and 
thus when we consider the general dynamics we cannot assume that T,ap = diag(Si, S2, S3) 
and Ra — 0. For our purposes, however, it suffices to consider an asymptotic spatial frame 
choice. In this paper we will use an asymptotic Fermi frame, i.e., but note that 

we cannot e.g. diagonalize the shear in such a frame, except asymptotically. On the Kasner 
subset /C, and the Jacobs subset JT", a Fermi frame choice leads to that the evolution 
equation (|16aP immediately yields 

Sa;9 = Sa^, (44) 

since g = 2 in both cases. Thus, because of that = = 1 — figtiff there exists a ellipsoidal 
ball (ellipsoid) of fix points in the (shear projected) J' (K.) case, which corresponds to a center 
manifold. However, in this case a temporally constant rotation of axes that diagonalizes E^^ 
so that Ea^ = diag(Ei, Ei, Ei), which leads to (|43)) . 

From ([43)) it follows that Aa, N^f), when a 7^ /?, and f2(i) {{i) 7^ stiff) are stable 
toward the past everywhere on and /C, with the exception of the non-transversally- 
hyperbolic so-called Taub points on /C, where (pi,P2,P3) = (1,0,0), and cycle (or equiv- 
alently (Ei, E2, E3) = (2, —1, —1), and cycle). However, we have shown that if lemma [3T2] 
holds, which depends on that the condition (IT^ holds, then A^ and both tend to zero 
toward the singularity and thus the states A^ — and — are past stable, even though 
they are not linearly stable everywherelll The decoupling of the equations (j43d|) from 
each other, and their shared linear stability properties in conjunction with the previous non- 
linear stability result, gives some support for the 'matter does not matter' conjecture (2), in 
the context of that the condition (IT^ holds, cf. also 151 [71 [T5]. 

The stability toward the past of Na depends on the sign of Pq = ^(1 + Eq). In the stiff 
case (ii) it follows from P5)) that the part of J that obeys Eq ~ — IVa, i.e. with E^, > —1 
or, equivalently, pa > 0, is stable w.r.t. Na perturbations toward the past; we will denote 
this part of as J''^. Outside J'^, J have an unstable mode associated with Na when 
Eq < —1, or, equivalently, when < 0. This follows from that only one of pi,P2, and p^ is 
negative, because pi + P2 + P3 — ^ and Pi + P2 + P3 = 1 ^ f^^stiff yields 

i(l - 2S) < p„ < i(l - E) < < i(l + t)<p^< i(l + 2E), (45) 

where E = \/l — ilstiff , and where (aPj) = (123), and cycle. Since we showed in Proposition 
13.31 that the past attractor in the stiff fluid case must be confined to we immediately get 
from requiring consistency with the stability analysis that it must be contained in the closure 
of the stable part of J^, i.e J'^. 

In the soft case (iii), it follows from (|43|) that K, is unstable everywhere toward the past 
(except at the excluded points Pmax = 1) with an unstable 7Vo,-mode when E^ < —1, or, 
equivalently, when < 0, since pi -I- P2 + Pa = 1 and Pi + P2 + P3 = 1 yields 

- I < < < p^ < I < p^ < 1, (46) 

where (a/Sj) = (123), and cycle. 

^^There exists other interesting choices, e.g. Ra = Eq S^-y, where (a/37) = (123), or cycle, and where 
Ea is equal to ±1; ea = (~lili~l) is connected with the Iwasawa frame used in e.g. |15| . ea = (1,1,1) 
is the frame choice used in |14l . For these choices Ra destabilizes parts of the fix point sets on J7, IC by 
inducing so-called frame transitions (also known as centrifugal bounces in a Hamiltonian context, see |15j). 
trajectories that connect one fix point representation of a type I solution with another, by means of an axes 
permutation |15] . 

There hence exists an intricate connection between avoidance of the Taub points, which in turn are part 
of the Taub subset described in Appendix [Bj and asymptotic locality via lemma 13.21 and condition I I12II . 
Determining exactly what this connection is poses a formidable and important challenge. In this context 
it is worth mentioning that there may exist an open set of solutions with so-called weak null singularities, 
which are not asymptotically silent or local; moreover, these singularities seem to be intimately associated 
with the Taub subset I18| . 
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The past instabilities in the stiff and soft cases are associated with so-called silent Bianchi 
type II curvature transitions, to use the nomenclature of [15] . i.e., orbits associated with 
Bianchi type II. We therefore take a closer look at this silent subset for the stiff and soft 
cases in a Fermi-propagated simultaneously diagonalized and Najs frame. 



3.5 The type II subset 

In the past asymptotic limit v^^^g = in the stiff case, and — in the stiff and soft cases, 
where i refers to a fluid with an asymptotically soft equation of state. We are thus interested 
in the subset on Bianchi type II that is described by a single stiff fluid with w^jfj = in 
the stiff case (ii), and the vacuum type II subset in the soft case (iii). We choose a Fermi- 
propagated shear eigenframe with = diag(Si, S2, S3), and project the dynamics onto 
Sct-fistiff-space, i.e., we disregard the test fields v^^y, in addition we set N^p = 0, except for 
a single component N^^ = N^, which we determine via the Gauss constraint, which yields 
= 12(1 - I]2 - rJstiff)- Eqs. (fT6al) and ((T6dl) then yield 

ao(2-S„) = -(2-q)(2-S„), (47a) 

ao(2-E^) = -(2-q)(2-I]/,), (47b) 

ao(4 + I]^) = -(2-g)(4 + I]^), (47c) 

^oatiff = -2(2 - <z)atiff , (47d) 

where (a/37) = (123), and cycle, and where q = 2(S^ + ilstiff), where fistiff = in the 
vacuum case. 

It follows that the solutions to (j47|) are trajectories that are straight lines when pro- 
jected onto Ea-space. Since — 2 < < 2, a = 1,2,3, and g < 2 on the type II subset, 
equations (j47|) show that E-y (Eq,,E^) is monotonically increasing (decreasing) toward the 
past and approaches a limit value on the type I boundary where q — 2. Together with 
the previous stability analysis, this implies that the solutions originate from Jacobi/Kasner 
fixed points with E^ < — 1 and end at Jacobi/Kasner fixed points with E^ > —1, when the 
direction of time is taken to be toward the singularity. Hence the global future attractor 
of (H71) is given by the fix points on /CU J" (/C in the vacuum case) for which E^ < —1, while 
the past attractor of (|T7)) is given by the fixed points on /C U J' (/C in the vacuum case) with 
E^ > —1. Using the nomenclature of [7j[T5], the solution trajectories are denoted as single 
curvature transitions, and they reflect and describe the outcome of the past instabilities 
associated with N^p described in the previous subsection. 



3.6 Past attractors on the silent boundary 

Combining the previous linear and non-linear stability analysis with the results in subsec- 
tion [331 notably the past asymptotic consequences for i>"-^ that were obtained via Eq. (|42|) . 
now allows us to make some firm statements about the past attractors on the silent bound- 
ary in cases (i) and (ii), and we also make some predictions about the past attractor in case 
(iii). Throughout we use an asymptotic Fermi frame. 

The ultra-stiff fluid case (i): 

Proposition 3.5. The ultra- stiff fluid case (i): If the past attractor A~]^^^^_^^-^^ is contained 
on the silent boundary, then it is given by 

■^ultra-stiff' ~ ' 

where is characterized by 

(E„^,7V„^,A„) = (0,0,0), aot = i, Q?ot = o, n«t = o, 

and 

^tot = ^^ultra-stiff = 1; <ltra-stiff = ^(i) = ^: 

^'W = (c2)(,) > i; v^.^ = uc^^) when {c^^)^i^ = i; v"^^^ = c^^^ when (c2)(,) < i. 
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Proof. This follows directly from Propositions 13.21 and I3.4[ and the results that followed 
from Eq. gl]). □ 



The stiff" fluid case (ii): 

Proposition 3.6. The stijf fluid case (ii): If the past attractor A~^^g is contained on the 
silent boundary, then it is given by 

A«ff = (J^r, (48) 

where (J^^)^ is characterized by 

Nats = 0, Aa= 0, (49a) 
— Sq,^, such that > ~1 {or, equivalently, pa > 0) Va, (49b) 

f^tot = atiff, Qtot^o, n^'^, = 0, (49c) 

where (I]i,S2,S3) = (3pi — l,3p2 — IjSps — 1) are the (non-ordered) eigenvalues ofHajS, 
and (apart from the stiff fluid(s) for which rJstiff — f^stiff = ^tot)- 

= 0; = 0; (50a) 

w^j) = when (c^)(j) > Pmax = 3(1 + Smax), (50b) 

^'^) = ^{i) (c^)(j) = Pinax = -1(1 + Sinax), (50c) 

^{i} = ^ti) '^^^'^ (Cs)(») < Pmax = 5(1 + Sniax), (50d) 

Proof. Follows directly from Proposition [231 Eq. (H^ . and the linear analysis of (|43p . Note 
that the velocity directions refer to the 'dominant' shear eigen-directions, as discussed after 
Eq. (021). □ 



The soft fluid case (iii): In the soft fluid case we have no proof, but the previous analysis 
suggests that there exists a past attractor A~^f^ subset, on the vacuum part of the silent 
boundary, that describes the asymptotic dynamics of a timeline in terms of Mixmaster like 
behavior (see [T2] sec. 6.4), where the asymptotic dynamics is approximated by an infinite 
heteroclinic sequence that reside on A~^f^ , which hence breaks asymptotic self-similarity ^21] ; 
A"oft is given by 

A-ft=^USr~, (51) 
where S^^ on the Kasner subset /C is described by 

— ^al3, (52) 

since we use a Fermi frame, and where g™cuum ^^iq silent vacuum Bianchi type II subset 

Since the conjectured attractor consists of the vacuum type I and II subsets, we expect 
that we (at least generically, recall footnote 8) can asymptotically diagonalize Eq,^ and A^q,^ 
(although the diagonalized shear directions will typically be different for different timelines). 
To describe the asymptotic dynamics we hence perform a constant rotation that diagonalizes 
an 'initial asymptotic' Kasner point SJ^^, which leads to that the subsequent dynamics is 
described by the same heteroclinic sequence (Kasner states joined by type II curvature 
transitions) as in the vacuum Bianchi type VIII and IX cases, see e.g. [12) . 

In the soft fluid case the description of -4~f^ also involves the asymptotic test fields 
v"-y Just as Eq,/3 and Nap oscillate perpetually, so do the fields v"-y The effects of a 
sequence of Kasner transitions by means of curvature type II transitions is two-fold: (a) a 
change of Kasner state, (b) a change of ordered Kasner shear eigen-directions associated with 
Pa < P/3 < P-y, where (a/37) — (123), or a permutation thereof. The curvature transitions 

^^This subset consists of the union of six disjoint Bianchi type II subset representations, each characterized 
by the sign of a single non-zero eigenvalue of Af^^^. 
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induce a sequence of 'tilt' (peculiar velocity) transitions. As the dynamics approach the 
attractor it follows that the asymptotic vacuum dynamics spend an increasing time near 
the Kasner fix points. This leads to that the asymptotic test fields have increasingly 
long periods of time to reach their past asymptotic states on the Kasner subset. This 
in turn implies that velocities for increasingly long times are either almost aligned, anti- 
aligned, or one or two of the associated speeds are zero, depending on (Cs)(i), (c^)^) {i ^ j), 
and Pmax, see subsection 13.31 However, this correlation is temporally broken and changed 
whenever there is a curvature transition. But since the curvature transitions are increasingly 
dominated in time by the Kasner states it follows that the probability of finding v"-^ and 
in the previously described correlated state increases with timeP^ 

4 Past stability and instability in the full state space 

We now turn to the discussion of the role of the past attractors on the silent boundary of 
the stiff and soft cases in the full physical state space. We have previously assumed that 
the dynamics approach the silent boundary toward the past and that i?Q* — > 0, and we have 
subsequently worked out the consequences of these assumptions. To check the consistency 
of this it is of interest to compute Ea^ 'on' the past attractors by inserting the attractor 
subset variable values in Fa^ in the evolutions equation d^Ea = Fa^Ep^, thus yielding 
a lowest order past attractor perturbation of E^ in the full state space. Since the past 
attractor resides on the type I subset in the stiff case (ii) and since we expect the type I 
subset to 'dominate' the Mixmaster dynamics in the soft case (iii) (the ultra-stiff case (i) 
has already been discussed previously), we insert the diagonalized shear values on Sf,lC in 
Fa^\ this yields the following equation for the individual Ea components: 

{E^T^doE^'\j^^ = 2-t^=^{l-p^), (53) 

and thus we see that E^ is stable toward the past everywhere on J , IC, except at the 
Taub points, as is to be expected, but which nevertheless yields support for the assumption 
Ej 0. 

Next we discuss Tq,. Recall that Tq, = —Ea^di \nH, and since -Eq.' — J> 0, then — > if 
we have chosen a gauge so that di IniJ does not blow up too fast. A way at looking at the 
evolution of Tq, is to heuristically regard the evolution equation (j78ap for asymptotically as 
an equation of the form dor a = cia^rfj+ba, where Oq,^ is F^^ computed on the past attractor, 
while ba is {da + Ua){<l + 1) calculated 'on' the past attractor, where E^ in da is computed 
by inserting the attractor values in F^^, which leads to the evolution equation d^Ea^ = 
o-a^ Ea^- This leads to that one can regard aa^ and ba as time dependent coefficients on 
a given timeline, effectively leading to an ODE for Tq, where the solution for is given 
by the solution to the homogenous equation d^ = aa^rp added to a particular solution 
associated with ba- However, due to that the homogenoeus equation for Tq, has the same 
character as that for it follows that the homogeneous solution tends to zero. Thus 
we require a gauge that is such that the particular solution also tends to zero, where the 
freedom in the gauge choice is reflected in the term 6^; we expect that we require a gauge 
such that [da + Ua){q + 1) tends to zero reasonably fast. Considering that g is 2 in the stiff 
case or 'almost always' 2 in the soft case due to 'Kasner dominance' this suggest that this 
is presumably a rather wide classlll For an example of a gauge with {Ma, Wa) = (0, 0) for 
which there is numerical support that — > (as well as Ua — >■ 0), see [3]. 

We now turn from considering reference congruences in general, to the issue if there 
are fluid congruences for which the Hubble normalized vorticity W" and acceleration U"' 

In |15| it was shown that one can expect that oscillations yield cumulative trends over time. This also 
pertains to the peculiar velocities, and the results in 1151 suggest that such trends depend on the stiffness of 
the equation of state; this is also suggested by numerical experiments for special models 1311 . 

'^^In \3\, where Ma = Wa = 0, it was noted that Ea^ = 0, Ua = yields an invariant boundary subset, 
where Va ^ leads to the same equations as those for spatially self-similar models. In [3], this subset was 
referred to as the silent boundary, but since flo r^ = (gSa^ — ^a^)ra r^ on this subset, which leads to that 
Ta — >■ 0, we have chosen to focus on the subset with ra = 0, which we here has referred to as the silent 
boundary. Furthermore, note that ra is stable toward the past on the 'extended' silent boundary. 
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will vanish asymptotically. Choosing the timelike reference congruence as one of the fluid 
congruences implies that for that fluid Va — 0, p = p, p — and Qa = = 0, while P — 
(again we drop the index (i)). The fluid equations reduce to (obtained by specializing 
the total matter equations ([771) Appendix |A] to a single comoving perfect fluid) 

don=[2q-l- 3w] n, (54a) 

O^cl (a„ - 2r„) O + (1 + w){Ua + ra)n, (54b) 

or equivalently, 

dQp= -3{p + p), (55a) 

= a„p+([7„+r„)(p + p). (55b) 

Assuming that the weak energy condition holds strictly for the fluid component at hand, 
i.e., p > and p + p > 0, makes it possible to introduce the particle density n and the 
chemical potential p, 

dn dp p+p dp dp 

— = — — . M= > — = — — . 56) 

n p + p n p p + p 

which, together with (|55p . yields 

do n ~ — 3n, (57a) 
0^ {da + Ua + rc)p, (57b) 

where a suitable function of n may be useful as a matter variable in the case w ^ const, 
see [22]. 

By applying do to (|57bp and using (|78al) and (|68aP we obtain 

do f>a - [Fj + (3c2 - 1 - q)Sj] tlfj + a„(3c2 - q). (58) 

Equations (1571) and (|78b[) together with applying (I68bp to In/i, and using the relation 
dlnp/dlnn — — dp /dp, yield 

^CjUp = {2,cl-q-l)W^, (59) 

which allows equation (|73bp to be written on the form 

do = {fJ + (3c2 - 1) 5 J + Wp. (60) 

Following Taub |23l[25, we let 

M=^, (61) 
where Mq = Mo{x^), which, via ([57]), ([75a|) . and ([75cl) yields that 

M, = M,ix')^m, (62) 
which gives that the time dependence of Ada is determined by i?a* since 

X„ = S„*M,. (63) 

Applying equations (|75ap and (|75cp to this result gives Wa — \M.Ej3''Cj^ Mi (a relation 
that is equivalent to the non- normalized coordinate frame expression = M^]^^Mj^/). Since 
Eq. ([S5)) implies that if E^ — > then M.a — > it remains to investigate if Wa and tends 
to zero toward the past. 

In the ultra-stiff case (i) it follows straight forwardly that Wa — > for the ultra-stiff 
fluid in the neighborhood of T. If in addition 5Q(3(Cs)uitra-stiff — ?) — > sufRciently fast, 
which can be shown to be a consistent condition by means of an analysis similar to that 
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of other isotropic singularities undertaken in [25.^ (see also [H]), then also Ua 0; this is 
to be expected since I'uitra-stiff' = asymptotically when measured some congruence that is 
assumed to satisfy the asymptotic surface formation condition (15a)) (we also expect that soft 
fluids with Cj > ^ in the ultra-stiff case satisfy the asymptotic surface formation condition 
since they lead to 0, cf. subsection 13.31) . 

Let us turn to the stiff (ii) and soft cases (iii). In analogy with subsection 13.41 let us 
study the stability of and f/^ by making a perturbation of J^/K. in a shear diagonalized 
Fermi frame. Eq. (j60p then yields 

W-^9o W^\j^^^ = 1 + 3c2 + E„ = 3(c2 (64) 

which requires + > Va in order for Wa — s- 0. 

On J'^, the stable satisfies > —1, Va {pa > 0, Va), and on this part Wa 
when Cg > 0. In the soft case (iii) mm{pi,p2,P3) — on K, and thus > ^ leads to that 
Wa everywhere on /C, but for fluids with < |, parts of the /C become past unstable 
with respect to the vorticity, and for dust {w = = 0) all of /C is unstable (except at the 
non-transversally-hyperbolic Taub points); hence the vorticity of dust does not vanish in 
the approach to the singularity. For < < ^ it is the cumulative effect over time of the 
factor + Pa that matters; to determine this effect would require a study by means of, for 
example, methods used in ;15^, which we will refrain from since it is not enough that the 
vorticity tends to zero in order for the asymptotic surface formation condition (|8ap to be 
fulfilled, it is also required that — !> 0. 

The analysis of (|58)) of the past asymptotic behavior of Ua is complicated by the term 
5a(3Cg— g). However, by considering its asymptotic expression, by inserting the asymptotics 
for q and c^, and by solving the evolution equation for 'on' the silent boundary (i.e., by 
perturbing the past attractor to lowest order), this term can be regarded as a time-dependent 
inhomogeneous term; similarly one can compute the factor before Ua on the r.h.s., which 
yields an equation of the form do Ua = aUa + where a and ba can be regarded as given 
time dependent functions on a given timeline. Hence the general solution can be obtained by 
adding a particular solution to the general solution of the homogeneous part, Ua — aUa- 
In order for the fluid to be asymptotically surface forming ([5a|) it is required that Ua 
generically, and a necessary condition for this is that [/q, — )> according to the homogeneous 
equation, which, when computed in a Fermi frame on jK,^ yields 

U-^doUa\j^X = - 1 - S„ = 3(c2 -Pa). (65) 

In the stiff case (ii) Ua ^ ^ requires that > Pmax on ^ i.e., the same condition as 
required for (which of course is to be expected). The condition that — s> holds 

everywhere on J'^, and in this sense holds for a generic solution, requires that — Pmax > 
everywhere on J''^, which leads to that cl — 1, i.e., none of the softer fluids will in this case 
fulfill condition ([5a|) . 

In the soft case (iii) there will always be parts of /C that are past unstable with respect to 
the fluid accelerations, the asymptotic behavior then depends on the cumulative effect of the 
factor (c^ —Pa), as it does with the vorticity; for fluids with < | all of K. is unstable in at 
least one mode, since one of the Pa > |; hence the acceleration does not vanish for any fluid 
with low sound speed. This in turn probably leads to that the condition on the particular 
solution for Tq, breaks down and hence Tq, does not tend to zero either. Hence comoving 
gauges for fluids with < | are not gauges that are compatible with the asymptotically 
surface forming condition (|8al) . and it may be that this is also the case when > | (to 
provide plausible arguments for this would require an extensive study using, for example, 
methods described in [H]). 

We conclude this section with a discussion of the asymptotic behavior of and p{i). 
The evolution equation for In on the silent class A subset is governed by the sign of the 
factor -[3 + vf.^ + E„/3 cf^^ cj^ vf^^] = -[2 + (1 - vf^^) + {2Sa0 + Sa0)c^,) cj) vf^^], see Eq. 
Because of Corollary p.ip . w^^j < 1 leads to that the factor associated with p(j-) is strictly 
negative, and hence it follows that lmi^o_^_^ p^j^ — >• oo, V i if a solution approaches the 
attractor on the silent boundary. 
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The evolution equation for In p^j) on the class A boundary is governed by the sign of the 
factor -[3 - v^^ - c^^^ cj^ vf^^] = -[3(1 - vf^^) + (25„^ - ^c.p)cf,^ cj) v^^], see Eq. m- 
Lemma [3.21 and w^^^ < 1 suggests that lima,o_j._oo — >■ C!0, V i, but unfortunately, we do 
not have a strict inequality in this case, and since needs to be evaluated in the interior 
physical state space it is not certain that the quantities that have been neglected in the above 
equation do not prevent — > oo; nevertheless oo is normally what is assumed in a 

BKL context and we do so here as well. Thus e.g. when appears in an equation that is 
used in an asymptotic context it refers to the limit when oo, and either represents 

an actual limit or a bound. 

5 Conclusions 

We have studied the past asymptotic dynamics of spacetimes with an arbitrary number of 
perfect fluids with non-zero peculiar velocities, and with general barotropic equations of 
state, where it has been assumed that the Hubble-normalized interactions can be asymptot- 
ically neglected. Using dynamical systems methods on a system of equations obtained from 
the 1+3 Hubble-normalized conformal orthonormal frame approach, we have reformulated 
two well known conjectures by Belinksii, Khalatnikov and Lifshitz, about properties at the 
vicinity of a generic spacelike cosmological singularity, to conditions on our variables ('the 
locality conjecture' (1) and 'the matter does not matter' conjecture (2)) and worked out 
the consequences of these assumptions. We have shown that from the assumption of 'the 
locality conjecture' (1) alone follows: 

• In the case where there exists at least one fluid with an equation of state that is ultra- 
stiff asymptotically to the past (i.e., the speed of sound Cs satisfies the inequality > 
1), then the Hubble-normalized shear and spatial curvature will vanish asymptotically 
along with the Hubble-normalized energy densities f2(i) of all the fluids but the one 
with the asymptotically stiffest equation of state, which will have a Hubble-normalized 
density parameter rJuitra-stiff of unity and a vanishing peculiar velocity. The peculiar 
velocities v"^-^ for fluids with rij-j) — asymptotically vanishes in this case when (c^)(j) < 
1/3 while = 1 if (c2)(,) > 1/3. 

• In the case where no fluid is asymptotically ultra-stiff, but at least one fluid is asymp- 
totically stiff {w — = 1), the past asymptotic temporal behavior is given by a 
Jacobs solution. In this case w^jfj = and ri(i) = asymptotically, Vi ^ stiff such 
that (Cs)(i) < 1. The pecuHar velocities v"--^ tend to zero when {c^)(i) > Pmax, 
where Pmax is the maximal shape parameter, while tj^^j — 1 when < Pmax. 

• In the case when there are no stiff or ultra-stiff fluids, but an arbitrary number of fluids 
with soft equations of state (cj. < 1), the past asymptotic state resides on the union 
of the Bianchi type I, II, VIq or VHq subsets on the silent boundary. Furthermore, in 
this case we have made the additional assumption that 'the matter does not matter' 
conjecture (2) holds, and have provided some arguments that the past attractor is 
Mixmaster like, with oscillations between Bianchi type I and II vacuum solutions on the 
silent boundary. The peculiar velocities of the fluids — which become test fleld — become 
forever oscillating, both in direction and amplitude as the oscillations change the shear 
via the so-called Kasner map (see e.g. [7]). Interestingly the extreme properties of 
the asymptotic spacetime geometry induces 'correlation effects' among the different 
peculiar velocities. 

By studying the vorticity and acceleration in a fluid comoving frame we come to the 
conclusion that the fluid comoving gauges for stiff and ultra-stiff fluids obey the 'locality 
conjecture' (1) and are therefore acceptable gauges. However, this is not the case for models 
with only fluids with asymptotically soft equations of state. In such models fluid comoving 
gauges for fluids with < 2/3, which notably include dust and radiation equations of state, 
are not compatible with the gauge requirement associated with 'the locality conjecture' (1), 
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and fluid comoving gauges for fluids with 2/3 < cj. < I may be inadmissible as well, although 
the latter is an open issue. 

In all, our results for spacetimes with multiple fluids agree with previous studies of special 
models with soft (e.g. [H [Ml EZ] ) , stiff (e.g. [2^), and ultra-stiff single fluid models [TT] . 
However, we here studied asymptotic dynamics in a general infinite-dimensional dynamical 
systems setting, where we pursued the consequences of BKL-like assumptions, and this led 
e.g. to the conclusion that comoving fluid gauges are, for the most physically interesting 
fluid cases, incompatible with BKL-like behavior. It follows that a matter element will 
always move w.r.t. to a frame that obeys the gauge requirements of 'the locality conjecture' 
(1), furthermore, even in the fluid comoving gauge a matter element will accelerate and 
pick up momentum w.r.t. the rest frame of the fluid. In this sense matter momentum 
will matter toward the singularity (which is not the case for fluids with stiff or ultra-stiff 
equations of state), even though 'matter does not matter' for the asymptotic spacetime 
geometry, answering a speculation posed in [M]. It was also beneficial to consider multiple 
fluids, since this made it possible to investigate the relative evolution of the fiuid themselves 
with some interesting results, like the dominance of the stiffest fluid, Eq. p2|) . and the 
peculiar velocity shear alignment in type I (Corollarv 13. 3p . which led to suggestive results 
about asymptotic correlations between different peculiar velocities for asymptotic oscillatory 
behavior, as discussed in subsection 13.61 

Our analysis has rested on the assumption that the silent boundary is approached, and 
in the soft fluid case on the further assumption of asymptotic vacuum dominance, and even 
though we found support for our assumptions it would be desirable for further study to 
establish firm results on all points. We therefore list three open problems. 

• There seems to exist an intricate connection between asymptotically approaching the 
Taub subset, described in Appendix [Bl and the violation of lemma and the con- 
dition (jl2p . Furthermore, in [15^ it was shown that one statistically with increasingly 
probability find the state of a solution in a small neighborhood of the Taub subset in 
the approach to an oscillating singularity. Moreover, there seems to be a connection 
between weak null singularities and the Taub subset [18]. Hence there is a need for 
a detailed separate study of the Taub subset, but such an analysis is unfortunately 
likely to pose a major challenge. 

• Another less formidable future possibility is to apply the methods in [15] to study 
cumulative trends for peculiar velocities, discussed in subsection 13.61 

• A third possible investigation would be to investigate what BKL-like conjectures, anal- 
ogous to the presently formulated ones, would imply for other sources. When does 
'matter does not matter' toward the initial singularity hold in this more general con- 
text? 
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A The 1+3 conformally Hubble- normalized dynamical 
systems approach 

To establish conventions and notation, we in this Appendix briefiy introduce the conformal 
1-1-3 Hubble-normalized dynamical systems approach. This constitutes a specialization of 
the results in [I], in combination with that we derive the general perfect fluid equations; for 
further details and motivation we refer to [T]. 
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In the conformal Hubble-normalized orthonormal frame approach, of. [TJ 129) , we intro- 
duce a conformal 'Hubble-normalized' orthonormal frame of g (or, equivalently, an orthonor- 
mal frame of G) according to g = H~^ G = H~^ rjab , where the one-forms Jl" are 
related to the conformal orthonormal vector fields da via ( $1", db) — S'^b- We align Bq with 
a timelike reference congruence, which leads to that do and da are given by: 

do = H-^ ea = M-^d^o, d^ = H'^ ^ Ma Mdo + E^'d^, (66) 

where A4 and Aia are the conformally Hubble-normalized threading lapse function and shift 
vector, respectively; denotes the time coordinate along the timelike reference congruence, 
while di = d^i, where are spatial coordinates {i = 1,2,3). Throughout we express 
the derivatives in all equations by means of the derivative operators do and da- Partial 
derivatives are thus 'weighted' with conformally normalized frame variables, and this is one 
of the main advantages of the present formalism. 

The deceleration parameter q and are objects that are kinematically defined by 

doH = -{q + l)H, daH = -raH. (67) 

For dimensional reasons, the above equations for the dimensional Hubble variable H, asso- 
ciated with the timelike reference congruence in the physical spacetime associated with g, 
must decouple from all equations that only involve dimensionless variables and operators. 

It is useful to write the dimensionless commutator equations on the following operator 
form: 

0^{da + Ua)do " {5 J do - fJ) dp, (68a) 
Q^2Wado-Cj dp, (68b) 

where 

fJ = -[n 5 J + + ej^ {W + R-*)] =qSj - - ej^ {W^ + i?^), (69a) 
Cj =ea^f^{d^~A^)~Nj, (69b) 

where we have used doH = — {q + 1) H to obtain the relationship q = —% = ^ I©, which 
relates the deceleration parameter q to the (Hubble-) conformal Hubble scalar H and ex- 
pansion 0; for the physical interpretation of the other quantities, see Section [1] 

We will be concerned with general relativity and hence we impose Einstein's field equa- 
tions: 

Gab = Tab, (70) 

where we have chosen c — 1 — SttG as units, where c is the speed of light in vacuum and G 
is Newton's gravitational constant. 

We make a 1-1-3 split of the stress-energy tensor Tab w.r.t. the tangential 4-velocity 
of the reference congruence in the physical spacetime g according to: 

Tab ^ pUaUb + 2q^aUh) + p hab + T^ab, (71a) 
hab = UaUb+ gab ^ hab u'' = 0; qa M° = 0, T^abu" = 0, 7r°a = 0, (71b) 

and hence the total stress-energy is encoded in the objects {p,p,qa,'i^ap)i where tt^q, = 0. 

The conformal transformation naturally yields new dimensionless matter variables by 
scaling p,p,qa,TTap with H^^, however, to conform with the standard definition il = 
p/(3i?^), we instead scale the matter variables as follows: 

{n, P, Q", Ua^} = {p, p, g", 7r„4/(3ff2). (72) 
The field equation for the dimensionless frame and commutator variables (obtained from 
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the commutator equations, the Jacobi identities, and the Einstein equations) are conve- 
niently grouped into evolution equations, and constraint equations: 

Evolution equations: 

doMa = fJ M0 + {da + Uc.)M-\ (73a) 

do = {fJ + qdj + 2^ J) Wp + ^Cj Up, (73b) 

doE^'^FjEp\ (73c) 

do ^ap = —(2 — 9)Sa,3 + 2e''"'(ct ^P)5 R-y - '^Sap + "S^afi - 2W^(q Rp) 

+ (a(„ + f7(„ + .!(„) ijp^ + 2(a(„ - r(„ + A(„) rp^ - Np^^ {Us + 2rs), 

(73d) 

do = Fj Ap + \{dp + Up){3qSj - Fj), (73e) 
do N"^ = (3(7<5^(" - 2F^(")iV^)^ + e^*("(a^ + U^)Fs^\ (73f) 

where 

i^a^ = g<5/-S/-e„^(VF^ + i?'^). (74) 

Constraint equations: 

= CjMp - 2X~iWa, (75a) 

= (a„ - f7„ - 2 Ac,) VK". (75b) 

= C„^£;^\ (75c) 

= 1 - I]2 - Ofe - O + iW^2 _ l^^^a _ i^2da - + r„) r", (75d) 

= (3,5„'' Ap + ej^ Nsp) - 3Q„ - (d^ + 2rp) - [C<,^ + 2e„'^^(f/^ + r^)] T^^ - 2r„, 

(75e) 

= iV^ - iS^ (e^/^^^^ + Nj) - {Fj - 2qSj + 2J:J) Wp, (75f) 
where = lE^^E"^, ^ p^^p^a^ and where 

Cj = ea'''^{d.y-A^)~Nj, (76a) 
g = 2E2 + 1 (f7 + 3P) „ 2 ^2 _ i^Q^ + U^~ 2(A„ - r„)] {U" + r"), (76b) 

- + 2e^^(„ iV^)5 A^ + d^{5\^ Ap^ + Np^s). (76c) 

r2fc = -^^; ^ = - 6^2 + 4da^", (76d) 
B^p = 2iV,^ 7V^^ - iVT^ A^„^, (76e) 

where A^ = The expression for q in (j76bl) was obtained from the Raychadhuri 

equation, which gives q its dynamical content; '^Sap, ^ can be interpreted as the trace- free 
and scalar parts, respectively, of the Hubble-normalized three-curvature, if the reference 
congruence is hypersurface forming (Wa = 0). The notation (...) stands for the trace-free 
part of a symmetric spatial tensor, i.e. ^(q/S) = A^p — ^Sap A^^. 

The equations (|73d|) . (|75dp . (|75ep . (|76bp . were all obtained from Einstein's field equations, 
and are thus dynamical in nature, furthermore, note that it is the total stress-energy content 
{ri, P, Q", Hap} that enters into these equations; all remaining equations were obtained 
from the commutator equations and the Jacobi identities, and are thus kinematical. If we 
want to stress that a quantity refers to the total stress-energy content below we will provide 
it with the subscript tot, e.g., fltot- 

These equations need to be supplemented with matter equations that depend on the 
chosen matter content, however, local conservation of the total energy-momentum yields 



-"^^In [l] there are sign errors in front of the terms ea^^ RpW-y, e"'^ (^^N^^jUs, and ea^^ RpW-y , in the 
equations that corresponds to ll75f|l . l)73d|l . and l)73b| l. respectively. 
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VfcT"'' = for the total ^"^ which for the 1+3 sphtted matter variables yieldsEZl 

Total matter equations: 

don = {2q-l)n-3P + 2A^ Q" - E^^n"^ - [da + 2([/„ + r„)] Q", (77a) 
do Qa = {fJ -{2-q) 6 J) Qp + (3(5„^ Ap + ej^ Nsp) n-^^ 

+ 2eJ^ Qfi - {df, + C//3 + 2rp) {P5j + IlJ) ~ Uo^n - r^^in - 3P). (77b) 

It is sonietimes useful to apply the commutator equations (|68p to log (77) and consider 
the following resulting auxiliary equations for Tq, (also possibly extending the above state 
space to include Va)- 

dorc = Fjr0 + {do, + Ua){q+l), (78a) 
= Cjrp^2{q+l)Wa. (78b) 

In this paper we assume that the matter consists of several perfect fluids with gen- 
eral barotropic equations of state. The stress-energy component of the z:th fluid satisfies, 
aT^i) = I^iy where I^^-^ represents the non-gravitational interaction term of the i:th fluid 
with the other fluids; since VqTjq^ = it follows that ^ ■ 1°-^ = 0. We are here going to 
assume that the Hubble-normalized interaction terms asymptotically tend to zero toward 
the singularity, and that the fluids, in this sense, are asymptotically non-interacting (this 
can still be the case even if the interaction energies tend to infinity). Using the Hubble- 
normalized version of the relation VaT^"^ = 0, since we assume that the Hubble-normalized 
interaction terms are asymptotically zero, leads to the following equations for D, and Va (to 
obtain less cumbersome expressions we drop the index (i)): 

dofl^ {2q~l-3w)n + [(3w - 1) - + 2(A„ -~ifc- r„) ~ 8^] Q" , (79a) 

dova - GZ^ [(1 - v^){3cl ~l~clA'^ vp) + (1 - cl){A^ + v^^) vp] 

- [Y.J + ej-^ {R^ + N^^ vs)] vp - A^v^ + tjm^ vp 

- {5 J - v^)Up - (1 -f w)-^{l - v2)[(l - w)5j - Awcl GZ^ vc,v^]rp 

^ [{Sj + 2cl GZ^ VaV^)d^{P5p'' + H^^) ~ (1 + cl)GZ^ dpQ^] , (79b) 

where G_ = 1 — c^v^, cj. — dp/ dp. In the above expressions all 'matter objects' refer to 
the i:th fluid component, except in q in (j79ap . since q obtains its dynamical content from 
the total source. A cosmological constant A can formally be regarded as a perfect fluid 
contribution with w = —1, which leads to the following Hubble-normalized stress-energy 
contribution: Ha — A/(3i?^) = —Pa, while = = H^^. Due to its definition and 
equation (|67|) . Ha satisfies daflA = 2(1 + q)flA, daH,A — 2raflA- 

B Invariant boundary subsets 

The physical interior of the state space is characterized by det(i?Q,') ^ 0, and, in the case 
of the interior of the perfect fluid state space, fitot 7^ 0. However, the asymptotes of most 
interior solutions reside on the boundaries of the interior state space, and hence it becomes 
necessary to study the dynamics on these boundaries as well. Some of these boundaries 
play a particularly important role. Notably we have the vacuum subset Qtot = (and hence 
^(i) — V i), and what we call the partially silent and the silent boundary subsets [T8] . 
The existence of these latter subsets is intimately connected with the homogeneity of (|73cp , 
which leads to the existence of boundary subsets of the interior subset (det(£'Q,*) ^ 0) such 
that the rank of the matrix E^^ is two, one, or zero. 




^^Notc that this is also a reasonable demand in the context of other metric theories than general relativity, 
i.e., H77| l has a broader area of application than the present general relativistic one. 
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Let us begin with the rank zero case. Our later discussion suggests that only a part of 
the subset = is of generic importance, namely the invariant boundary subset 

{E^\Ma.,W^,Uo.,r^)^0, (80) 

see ([73]) - ([77|) . which we denote as the silent boundary, where A^q, = 0, = yields da — 

0. On this subset, described by a state vector S (see Eq. ([T5|) for the case of several perfect 
fluids), there exists a coupled set of ordinary differential equations and algebraic constraints 
that are identical to those of spatially homogeneous models. This can be seen as follows. In 
the spatially homogeneous case a spatially homogeneous foliation with orthogonal timelines 
{Ma = Wa = 0) leads to (M,i?,S) = {M{x°), H{x°),S{x°)), and hence [>„ = r„ = and 
SqS = Ea^di S = 0, and as a consequence the equations for Ea^ {det{Ea^) ^ 0) decouple 
from the rest of the variables in S, and thus one often only considers the equations for the 
'essential' variables of the state vector S, cf. [30]. Although the equations for S coincide for 
the spatially homogeneous case and the silent boundary, there is a fundamental difference; 
in the spatially homogeneous case the constants of integration are really constants, but on 
the silent boundary the integration coefficients are spatial functions, since the state space 
in this case corresponds to an infinite set of identical copies — one for each spatial point. 

A similar phenomenon happens when the rank of the matrix Ea^ is one or two, which 
leads to boundary subsets on which the dynamics is identical to that of models with spatial 
symmetry orbits of dimensions two or one, respectively. We refer to these subsets as partially 
silent boundaries; in these cases there are two or one spatial coordinates, respectively, that 
act as an index set, in analogy with what happens for the state vector in the silent boundary 
case. 

Yet another, overlapping, boundary is of interest — the Minkowski subset. In the present 
formulation this subset corresponds to the Minkowski solution/spacetime in foliations for 
which H > 0. Hence it is characterized by that the Hubble-normalized curvature is zero, 

1. e., both r^tot and the Hubble-normalized Weyl tensor are zero. 

There are many subsets on the Minkowski boundary, but one seems to be of particular 
importance, the Taub subset (so denoted because it is related to the Taub representation 
of the Minkowski spacetime), which we define as a subset that, in addition to f^tot = 0, 
satisfies ^ = 0, i.e., V,k = 0, and {Ma,Wa,Ua,ra) = 0, and hence = 1 and q — 2. 
Furthermore, these conditions implies det(SQ^) = 2, and that it is possible to introduce a 
Fermi propagated frame in which 'Eap = diag(2, —1, —1), or cycle. 
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